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ABSTRACT 

A Presburger formula is a boolean formula with variables 
in N that can be written using addition, comparison (<, =, 
etc.), boolean operations (and, or, not), and quantifiers (V 
and 3). We characterize sets that can be defined by a Pres- 
burger formula as exactly the sets that can be represented by 
a rational generating function; a geometric characterization 
of such sets is also given. In addition, if p = (pi, . . . ,p„) are 
a subset of the unbound variables in a Presburger formula, 
we can define a counting function g{p) to be the number of 
solutions to the formula, for a given p. We show that every 
counting function obtained in this way may be represented 
as, equivalently, either a piecewise quasi-polynomial or a 
rational generating function. Finally, we translate known 
computational complexity results into this setting and dis- 
cuss open directions. 

Categories and Subject Descriptors 

F.4.3 [Mathematical Logic and Formal Languages]: 

Formal Languages — decision problems; G.2.1 [Discrete Math- 
ematics]: Combinatorics — counting problems, generating func- 
tions; F.2.2 [Analysis of Algorithms and Problem Com- 
plexity]: Nonnumerical Algorithms and Problems — Geo- 
metrical problems and computations 

General Terms 

Theory, Algorithms 

Keywords 

Discrete geometry, Ehrhart polynomials, generating func- 
tions, Presburger arithmetic, quasi-polynomials, rational gen- 
erating functions 

1. INTRODUCTION 

A broad and interesting class of sets are those that can be 
defined over N = {0, 1,2,.. .} with first order logic and ad- 
dition. 



Definition 1. A Presburger formula is a boolean formula 
with variables in N that can be written using addition, com- 
parison (<, =, etc.), boolean operations (and, or, not), and 
quantifiers (V and 3). We will denote a generic Presburger 
formula as _F(b;u), where b are the bound variables (those 
associated with a quantifier) and u are the unbound vari- 
ables; we use bold notation like b to indicate vectors of 
variables. 

We say that a set S C N"* is a Presburger set if there exists 
a Presburger formula F{h; u) such that 5 = {u € N"^ : 
nb;u)}. 



Example 1. The Presburger formula 

^^(6, It) = (u > 1 and 36 G N : b + b + l^u) 

defines the Presburger set {3, 5, 7, . . .}. Since multiplica- 
tion by an integer is the same as repeated addition, we 
can conceive of a Presburger formula as a boolean combi- 
nation of integral linear (in)equalities, appropriately quan- 
tified: 3b (u>l and 2b+l = u]. 



Presburger proved [33] that the truth of a Presburger sen- 
tence F{h) is decidable. In contrast, a broader class of sen- 
tences, where multiplication of variables is allowed, is unde- 
cidable; this is the negative solution to Hilbert's 10th prob- 
lem, given by Davis, Putnam, Robinson, and Matiyasevich 
(see, for example, |17|). 

We would like to understand more clearly the structure of a 
given Presburger set. One way to attempt to do this is to 
encode the elements of the set into a generating function. 



Definition 2. Given a set S C 
ing function is 



its associated generat- 



/(S;x) = ^x== Yl ' 
ses (si,...,si)ss 



For example, if S is the set defined by Example [T] then 
f{S; x) = x'' + x' + x' + ■■■ = 

1 — x^ 

We see that, in this instance, the generating function has a 
nice form; this is not a coincidence. 



Definition 3. A rational generating function is a function 
that can be written in the form 

g(x) 

(I-Xbl)--.(I-Xbfc)' 

where g(x) is a polynomial in Q[x] and G N'' \ {0}. 

We will prove that 5 C N'' is a Presburger set if and only if 
f{S; x) is a rational generating function. These are Proper- 
ties 1 and 3 in the following theorem: 

Theorem 1. Given a set S C N'', the following are equiv- 
alent: 

1. S is a Presburger set, 

2. S is a finite union of sets of the form Pn(A + A), where 
P is a polyhedron, A £ Z'*, and K^-U^ is a lattice. 

3. /(S';x) is a rational generating function. 

Property 2 gives a nice geometric characterization of Pres- 
burger sets; the set in Example [1] can be written as [3, oo) n 
(1 + 2Z). 

We are particularly interested in generating functions be- 
cause of their powerful flexibility: we can use algebraic ma- 
nipulations to answer questions about the set. For exam- 
ple, f{S; 1, 1, . . . , 1) is exactly the cardinality of S (if finite). 
More generally, we may want to count solutions to a Pres- 
burger formula as a function of several parameter variables: 

Definition 4- The Presburger counting function for a Pres- 
burger formula F(h; c, p) is 

5F(p) = #{ceN^: F(b;c,p)}. 

Note that c (the counted variables) and p (the parameter 
variables) are unbound variables. 

A classic example is to take F{c,p) to be the conjunction of 
linear inequalities of the form aici + ■ ■ ■ + OdCd < oop, where 
fli G Z. Then gpip) counts the number of integer points in 
the p^^ dilate of a polyhedron. 

Example 2. If F{ci, C2,p) is 2ci -I- 2c2 < p, then the set 
of solutions (ci,C2) £ lies in the triangle with vertices 
(0,0), (0,p/2), (p/2,0), and 

-'"4([IJ-')([|J+^) 




+ |p + 1 if p is even, 
+ |p + I if p is odd. 



The nice form of this function is also not a coincidence. For 
this particular type of Presburger formula, Ehrhart proved 
[19| that the counting functions are quasi-polynomials: 



Definition 5. A quasi-polynomial (over Q) is a function 
g : N" —J- Q such that there exists an n-dimensional lattice 
A C Z" together with polynomials gx(p) G Q[p\j one for 
each A G Z"/A, such that 

fl(p) = ^(P), for p G A. 

In Example O we can take the lattice A — 2Z and each 
coset (the evens and the odds) has its associated polynomial. 
We need something slightly more general to account for all 
Presburger counting functions: 

Definition 6. A piecewise quasi-polynomial is a function 
g : N" — !■ Q such that there exists a finite partition [J^{Pi Pi 
N") of N" with Pi polyhedra (which may not all be full- 
dimensional) and there exist quasi-polynomials gi such that 

9{p)=g^{p) forpGP, nN". 

One last thing that is not a coincidence: For the triangle in 
Example (2] we can compute 

'^^gF{p)x^ = 1 -\- X -\- 3x^ -\- 3x^ -f 6x'^ + ■ ■ ■ 

1 

" (l-x)(l-x2)2' 

a rational generating function! The following theorem says 
that these ideas are - almost - equivalent. 

Theorem 2. Given a function g : N" — ^ Q and the fol- 
lowing three possible properties: 

A. g is a Presburger counting function, 

B. g is a piecewise quasi-polynomial, and 

C. X^psN" gi'P)'^^ is 0, rational generating function, 

we have the implications 

Remark 1. Proving Theorem [2] will give us much of The- 
orem [T] using the following idea. A set S C Z'' corresponds 
exactly to its indicator function 

, , I 1 if u G 5, 
^^(")^\0 iiuiS. 

If S is a Presburger set defined by F{h; u), then 

Xs(u) = #{c G N : F(b; u) and c = 0} 

is a Presburger counting function. 

In light of Theorem [T] we might wonder if there is a sense 
in which B ^ A. Of course we would have to restrict g, 
for example requiring that its range be in N (Theorem [T] 
essentially restricts the range of g to {0, 1}, as it must be 



an indicator function). The implication still does not hold, 
however. For example, suppose the polynomial 

g{s,t) = {t- s'^f 

were a Presburger counting function given by a Presburger 
formula F{h; c, s, t), that is, 

g{s,t) = #{cen'' : F{h;c,s,t)}. 

Then the set 

{{s,t)eN^: $c F{h;c,s,t)} = {{s,t) en'' : g{s,t) ^ 0} 

= {{s,s^) : s en} 

would be a Presburger set. This is not the case, however, as 
it does not satisfy Property 2 in Theorem[l] If the parameter 
is univariate, however, the following proposition shows that 
we do have the implication B ^ A. 

Proposition 3. Given a function g -. N Q, if g is a 
piecewise quasi-polynomial whose range is in N, then g is a 
Presburger counting function. 

In Section 4, we prove Theorem [T] Theorem [5J and Propo- 
sition |3l Theorem [1] was originally proved in the author's 
thesis ^B]; in this paper, it is put into context as a con- 
sequence of the more general Theorem (2] In Section 2, we 
briefly survey related work. In Section 3, we present the pri- 
mary tools we need for the proofs. In Section 5, we turn to 
computational questions; this section surveys known results, 
but restates them in terms of Presburger arithmetic. 

2. RELATED WORK 

Presburger arithmetic is a classical first order theory of logic, 
proven decidable by Presburger [33]. Various upper and 
lower bounds on the complexity of decision algorithms for 
the general theory have occupied the theoretical computer 
science community, see 8, 16, 20, 22, 23, 32 . 

A finite automata approach to Presburger arithmetic was 
pioneered in [121 114| . and continues to be an active area of 
research (see, for example, [101 1151 \45\ ). This approach 
is quite different from the present paper's. 

The importance of understanding Presburger Arithmetic is 
highlighted by the fact that many problems in computer sci- 
ence and mathematics can be phrased in this language: for 
example, integer programming [29II38[ . geometry of numbers 
[131 127| , Grobner bases and algebraic integer programming 
[411 143[ , neighborhood complexes and test sets [35 , 42 , the 
Frobenius problem ,34,, Ehrhart theory [T] [19], monomial 
ideals [31J , and toric varieties [2T] . Several of the above refer- 
ences analyze the computational complexity of their specific 
problem. In most of the above references, the connection to 
Presburger arithmetic is only implicit. 

Connections between subclasses of Presburger arithmetic 
and generating functions are made explicit in [3] [4| [5] . Con- 
nections between rational generating functions and quasi- 
polynomials have been made in [191 1391 140] . and the al- 
gorithmic complexity of their relationship was examined in 
[44[ . This current paper is the first to state and prove a 



general connection between Presburger arithmetic, quasi- 
polynomials, and rational generating functions. 

The algorithmic complexity of specific rational generating 
function problems has been addressed in, for example, [T] 
[5l |9] I18l 1241 126] . Several of these results are summarized in 
Section 5 of this current paper. 

3. PRIMARY BACKGROUND THEOREMS 

Here we detail several tools we need for the proofs in Section 

4. When these tools have algorithmic (polynomial-time) ver- 
sions, we mention them too, which will help our discussion 
in Section 5. 

The first tool we need is a way to simplify Presburger formu- 
las. As originally proved [33] by Presburger (see [32] for a 
nice exposition) , we can completely eliminate the quantifiers 
if we are allowed to also use modular arithmetic. 

Definition 7. An extended Presburger formula is a boolean 
formula with variables in N that can be written using any- 
thing from Presburger arithmetic, in addition to statements 
of the form 

aiyi + 022/2 + ■ ■ ■ + aiyi = ao (mod m), 

where Oi e Z and m e Z+. 

Theorem 4. Given a formula F{h; u) in extended Pres- 
burger arithmetic ( and hence any formula m Presburger arith- 
metic), there exists an equivalent quantifier-free /orrottZa G(u) 
such that 

{u G : F(b; u)} = {u e : G(u)}. 

The set from Example [T] can be written without quantifiers 
as 

M > 1 and u = 1 (mod 2) . 

Next, we give two theorems that tie in generating functions. 
The first gives us a way to convert from a specific type of 
Presburger set to a generating function. 

Theorem 5. Given a point X e if , a lattice A C Z'*, and 
a rational polyhedron P C R>q, /(Pn(A-l-A); x) is a rational 
generating function. In addition, for fixed d, this rational 
generating function can be found in polynomial time. 

The first step to proving this is to use Brion's Theorem [11] . 
which says that the generating function can be decomposed 
into functions of the form f(^K n (A -|- A);x), where K is 
a cone. Then, notice that integer points in cones have a 
natural structure that can be encoded as geometric series: 

Example 3. Let iiT C be the cone with vertex at the 
origin and extreme rays u = (1,0) and v = (1,2). Using the 
fact that the lattice (uZ -I- uZ) has index 2 in Z^, with coset 
representatives (0, 0) and (1, 1), every integer point in K can 



be written as either (0, 0) + Aiu + A2V or (1, 1) + Aiu + A2V, 
where Ai, A2 G N. Therefore 

f{K n Z2;x) = (x*"'") + x(i''))(l + x" + x^" + ■ ■ • ) 

■ (1 + x" + x^" + • • • ) 

x(o,o)^x{i,i) 



(1 - x")(l - x^) 



See O Chapter VIII], for example, for more details. The 
complexity version of Theorem [S] is due to Barvinok [I]. For 
polynomial-time algorithms yielding generating functions, 
the output will generally be represented as a sum of the 
basic rational functions given in Definition [S] getting a com- 
mon denominator and simplifying into one fraction may take 
exponential time. 

Next, we would like to be able to perform substitutions on 
the variables in a rational generating function and still retain 
a rational generating function; particularly, we would like to 
substitute in I's for several of the variables. 



Theorem 6. Given a rational generating function /(x), 
then 

g(z) = /(z'\z'=,...,z''^), 

with li G N'^ , is also a rational generating function, assuming 
the substituted values do not lie entirely in the poles of f . 
In particular, substituting in Xi — z =1 yields a rational 
function, if 1 is not a pole of f. 

Furthermore, for a fixed bound on the number of binomials in 
the denominator of f, there is a polynomial time algorithm 
to compute g as a rational generating function. 



The existence version of this theorem is immediate: simply 
substitute in. The complexity version is more difficult and 
is due to Barvinok and Woods [5]. The problem is that we 
may have the rational generating function input as a sum 
of rational functions. Unfortunately, the substituted values 
may be in the poles of some of the terms of the sum without 
being in the poles of the entire generating function; careful 
limits must be taken. 



Example 4- Suppose we would like to substitute x = 1 
into 

1 -.1000 

f^^^-ih-h-x 

Combining and simplifying to f{x) = 1 -I- a; -I- • • • -I- x^^^ is 
horribly inefficient. Since a; = 1 is a pole of both fractions, 
we must be careful: Consider the respective Laurent series 
expansions of the terms: 



: C-liX 



{x - 1)-' + coix - 1)° + ci{x~ + ■ 



l~x 

= d-i{x ~ 1) -|-do(a; — 1) + di{x - 



X 



+ ■ 



0. 



We know that 1 is not a pole of f{x), so c_i -I- d_i 
Substituting in a; = 1 into the sum, we will get simply co+do 
Carefully calculating co and do will give us /(I) = 1000. 



Finally, we need a connection between Presburger formulas 
and quasi-polynomials. This is given by Sturmfels [40] : 



Definition 8. Given ai, . . . , a^j G N", the vector partition 
function g : N" N is defined by 



5(p) = #{(Ai,...,Ad)G: 



p = Aiai H h Ada^}, 



that is, the number of ways to partition the vector p into 
parts taken from {a^}. 



Theorem 7. Any vector partition function is a piecewise 
quasi-polynomial. 



This is a surprising result. See [6] for a self-contained expla- 
nation utilizing the partial fraction expansion of the gener- 
ating function 



E 5(P) 

pgN" 



1 



(1-X«l)---(1 



this equality can be obtained by rewriting the rational func- 
tion as a product of infinite geometric series. For example, if 
ai = 1 and a2 — — 2, then the vector partition function 
is encoded by the generating function 

1 



(l-a;)(l-x2)2' 

We saw previously that this generating function corresponds 
to the quasi-polynomial in Example (2] 

A complexity version of Theorem [7| is discussed in Section 
15.21 an appropriate representation for piecewise quasi-poly- 
nomials must first be decided on. 



4. PROOFS 

4.1 Proof of Theorem H 



Given a Presburger counting function, g(p) = #{c G N'' : 
F{h; c, p)}, we first apply Presburger Elimination (Theorem 
lU to J" to obtain a quantifier free formula, G(c,p), in ex- 
tended Presburger arithmetic such that g{p) = #{c G N'' : 
G(c, p)}. Integers which satisfy a statement of the form 

aipi + ■ ■ ■ + OnPn + On+iCi + ■ ■ ■ + a„+dCd = flQ (mod m) 

are exactly sets A + A, where A G Z"^'' and A is a lattice 
in Z"'^''. Since G(c, p) is a Boolean combination of linear 
inequalities and these linear congruences, we may write the 
set, S, of points (c,p) which satisfy G(c, p) as a disjoint 
union 



where, for 1 < i < k, Pi C is a polyhedron, Ai is a 

sublattice of 2"+"^, and Ai is in Z"+'^. (To see this, convert 
the formula into disjunctive normal form; each conjunction 
will be of this form Pi fl (Ai -|- Ai); these sets may overlap, 
but their overlap will also be of this form.) 



Let Si = Pi n (Ai + Ai). By Theorem O we know we can 
write f{Si;y,x.) as a rational generating function, and so 

/(S;y,x) = ^/(S,;y,x)= ^ y^x^ 

i (c,p): G(c,p) 

can be written as a rational generating function. Finally, we 
substitute y = (1, 1, . . . , 1), using Theorem |6] to obtain the 
rational generating function 

^#{cGN^ G(c,p)}xP = X]g(p)xP. 
p p 



C ^ B. 

It suffices to prove this for functions g such that g{p):>i^ 
is a rational generating function of the form 



(1 - X=i)(l - X'^a) ■ • ■ (1 - X'^fc) ' 

where q G N",ai G N" \ {0}, because the property of be- 
ing a piecewise quasi-polynomial is preserved under linear 
combinations. Furthermore, we may take q = (0, 0, . . . , 0), 
because multiplying by x'^ only shifts the domain of the 
function g. Expanding this rational generating function as 
a product of infinite geometric series, 

J2 5(p)xP = (1 + ai + a? + • • ■ ) • ■ • (1 + afc + ai + ■ • • ), 
p 

and we see that 

g{p) = . . . , Afc) e N'= : p = Aiai + ■ ■ ■ + Afcafc}. 

This is exactly a vector partition function, which Theorem 
[7] tells us is a piecewise quasi-polynomial. 



B ^ C. 

Any piecewise quasi-polynomial can be written as a linear 
combination of functions of the form 

g(p)^/p" if peJ'nlA + A), 

1 otherwise, 

where a G N", P C R^g is a polyhedron, A G Z", and A 
is a sublattice of Z". Since linear combinations of ratio- 
nal generating functions are rational generating functions, 
it suffices to prove it for such a g. Let Cij, for 1 < i < n and 
1 < i < ii) be variables, and define the polyhedron 

Q = {(p,c) GN"+'"i + -+''" : 

p G P and 1 < Cij < pi for all Cij}. 

This Q is defined so that ^^{c : (p, c) G Q} is p^^ ■ ■ -p"" = 
p'* for p G P (and otherwise). Using Theorem [5] we can 
find the generating function for the set Q n (A + A) as a 
rational generating function. Substituting c = (1, 1, ... , 1), 
using Theorem[Sl gives us X]p5(p)^^ ^ * rational generat- 
ing function. 

4.2 Proof of Theorem [1] 

Given a set 5* C Z"^, define the indicator function, xs : — > 
{0, 1}, as in Remark[T] Define a new property: 



2'. xs is a piecewise quasi-polynomial. 

Translating Theorem [2] into properties of 5* and xs, we have 

1 ^ (2' « 3). 
So we need to prove 2 1 and 2' ^ 2. 

2^1. 

This is straightforward: the property of being an element of 
A -I- A can be written using linear congruences and existential 
quantifiers, and the property of being an element of P can 
be written as a set of linear inequalities. 



2' ^ 2. 

Since xs is a piecewise quasi-polynomial, it is constituted 
from associated polynomials. Let us examine such a poly- 
nomial g(p) that agrees with xs on some Pn(A-l-A), where 
P C is a polyhedron, A G Z", and A a sublattice of Z". 
It suffices to prove that 2 holds for S fl P fl (A -I- A), since S 
is the disjoint union of such pieces. 

Ideally, we would like to argue that, since q only takes on 
the values and 1, the polynomial q must be constant on 
P n (A -I- A), at least if P is unbounded. This is not quite 
true; for example, if 

P = {{x,y) eR^ : a; > and < 2/ < l}, 

then the polynomial q{x,y) = y is 1 for y = 1 and for 

What we can say is that q must be constant on any infinite 
ray contained in P n (A + A) : if we parametrize the ray by 
x(f) — {xi{t), ■ ■ ■ ,Xn{t)), then (j'(x(t)) is a univariate poly- 
nomial that is either or 1 at an infinite number of points, 
and so must be constant. Inductively, we can similarly show 
that q must be constant on any cone contained in P. 

Let K be the cone with vertex at the origin 

K = {y eR" : y + P c P}. 

Then K is the largest cone such that the cones ■k + K are 
contained in P, for all x G P; AT is often called the recession 
cone or characteristic cone of P (see Section 8.2 of [37)). 
and the polyhedron P can be decomposed into a Minkowski 
sum K + Q, where Q is a bounded polyhedron. We can write 
P n (A -I- A) as a finite union (possibly with overlap) of sets 
of the form 

Qj = {vj+K)n{\ + A), 

for some Vj, and on each of these pieces q must be constant. 
If q is the constant 1 on Qj, then Qj is contained in 5", and 
if q is the constant 0, then none of Qj is in S. Since S is a 
finite union of the appropriate Qj, S has the form needed 
for Property 2. 

4.3 Proof of Proposition |3] 

Given that g is a piecewise quasi-polynomial with range in 
N, we must find a Frobenius formula F{h;c,p) such that 



dip) = #{c G N"* : F(b; c,p)}. It suffices to find an F that 
agrees with g for sufficiently large p, because for any finite 
set {pi\ we may include 

(p = Pi) ^ (1 < ci < g{pi) A C2 = • ■ • = Cd = O) 

in the creation of F, so that the number of c satisfying F is 
exactly g{pi). 

Since the domain of g is one-dimensional, a piecewise quasi- 
polynomial can be thought of as a function that is a quasi- 
polynomial for sufficiently large p. Therefore we may as- 
sume, without loss of generality, that g is a quasi-polynomial. 
Each lattice coset may be handled separately and combined 
to form the final formula, F. Therefore we may assume that 
y is a polynomial. 

Let m be the least common multiple of the denominators 
of the coefficients of jr € Q[p]. Suppose p = mh -\- i, for 
some 6 (which can be encoded in a Presburger formula, with 
m and i constants and h a bound variable). Considering 
[mh + i)^ as a polynomial in h, all coefficients but (possibly) 
the constant term are multiples of m, so considering g{mb + 
i) as a polynomial in b, all coefficients but (possibly) the 
constant term will be integral. Since g{mb + i) is integer- 
valued, the constant term must also be an integer, that is, 
g{mb + i) £ Z[b]. Since each residue class of p mod m can 
be handled separately and combined at the end, we may 
assume without loss of generality that m = 1 and g £ Z[p]. 

Let g{p) = X]f=o ^ ^- Since g{p) > for all p, the 

leading coefficient, a^, must be positive. Then afep* is the 
number of c £ N'''+^ such that 

(1 < Co < afc) A (1 < ci < p) A ■ ■ ■ A (1 < Cfc < p). 

All other positive terms of g can be similarly encoded as 
counting solutions to a Presburger formula (carefully defined 
so the solution sets are disjoint). To take care of a negative 
term Oip*, take p sufficiently large so that a^p'' sufficiently 
dominates all other terms (as before, the finite set of smaller 
p can be dealt with separately); now remove a^p' points from 
the set of solutions. For example, 

^(^(co = Ci+2 = • • • = Cfc = 1) A (1 < ci < a,) 

A (1 < C2 < p) A ■ ■ ■ A (1 < c,+ i < p)) 

removes a^p' solutions, as long as p > a; (to do this for 
multiple i, define these subtractions carefully so that the 
solution sets do not overlap). All together, this give us our 
desired _F(b;c,p) certifying that g is a Presburger counting 
function. 

5. COMPUTATIONAL ASPECTS 

We first describe computational aspects oi A ^ {B and C). 
Many partial results are known here, but there are several 
open problems. Then we discuss B C, which works 
nicely: not only are B and C logically equivalent, but, af- 
ter we make sense of the problem, they are computationally 
equivalent (in fixed dimension). Finally, we discuss other 
related complexity results. 

5.1 {B and C) 



Note that another reasonable class of problems to study 
would be sentences where we are also allowed multiplication 
of variables, for example 

3a G N, 36 G N : + 2fe^ < 31. 

In general, however, these problems are very hard. In fact, 
there is a certain multivariate polynomial p{xo,xi, . . . ,Xd) 
such that the class of problems 

Given a G N, decide whether 3fei , 3fe2 , . ■ . , 36^ : 
p{a,bi,b2,...,bd) =0, with 6, G N 

is undecidable. This is a consequence of the DPRM-theorem 
(after Davis, Putnam, Robinson, and Matiyasevich, see, for 
example, [TT]), which solves Hilbert's 10th problem in the 
negative. Hilbert asked for an algorithm that, given a 
multivariate polynomial p, would decide whether p has any 
integer roots. 

The Presburger arithmetic problem (deciding whether a given 
F{h) is true or false) is, at least, decidable, as originally 
proved [33) by Presburger in 1929. Since then, better algo- 
rithms have been found. For example, D. Oppen gave an 
algorithm [32], based on work of D. Cooper [16) . with run- 

ning time 2 , where <j) is the input size of the problem and 
c is a constant. Nevertheless, Fischer and Rabin gave lower 
bounds on the running time [20) : any algorithm that solves 
all Presburger arithmetic problems will sometimes take at 

least 2^ * steps, where c' is a constant. 

Among the simplest Presburger sentences are integer pro- 
gramming problems, which are themselves NP-complete if 
the dimension is not fixed. This suggests that any interest- 
ing polynomial time results will require a fixed number of 
variables. 

We discuss three types of problems for a Presburger formula 
F(b;u): 

• The decision problem: Is 3u _F(b;u) true or false? 

• The counting problem: How many solutions u to F{h; u) 
are there? 

• The generating function prohlem: Compute ^(b.^j x" 
as a rational generating function. 

Clearly, the counting problem is harder than the decision 
problem. Now notice that if j^j^.^j x" can be computed 
in polynomial time, then the counting problem #u -F(b; u) 
can be solved in polynomial time: use Theorem [^ to plug 
in X = (1,1,... ,1). Consequently, the generating func- 
tion problem is the hardest of the three and the decision 
problem 3u F{h; u) is the easiest. Historically, polynomial 
time results have been discovered first for the decision prob- 
lems, and these solutions have not needed generating func- 
tions. However, polynomial time results for counting prob- 
lems have generally required the full power of generating 
functions. We now discuss several of these results. 

For a fixed number of variables, deciding 3u G(u) , where G 
is quantifier-free, can be accomplished in polynomial time. 
This is effectively Lenstra's algorithm for integer program- 
ming |29) , which decides whether there exists a u satisfying 



the conjunction of several linear inequalities. The only ad- 
ditional step we need is to prove that a boolean combination 
of linear inequalities can be put in disjunctive normal form 
in polynomial time. This is slightly surprising, because it 
is not doable for general boolean functions: for example, 
the disjunctive normal form for V ^,2) requires 2" 

disjunctions. 



Proposition 8. Fix d. There is a polynomial time algo- 
rithm which, given a quantifier-free formula G{ui ,U2,...,Ud) 
consisting of linear inequalities and boolean operations, con- 
verts G into disjunctive normal form. 

Proof. The inequalities appearing in G cut R'' into many 
polyhedral pieces. It might appear at first glance that the 
number of such pieces could be 2^, where is the total 
number of inequalities in G. Nevertheless, the number of 
pieces is bounded by 



$(d,N) 



+ 



N 



See Section 6.1 of [30] for a proof by induction. We have that 
$(d, N) is a polynomial in TV of degree at most d (where d 
is fixed), and following the inductive proof, we see that the 
description of each piece (and the lower-dimensional polyhe- 
dral pieces on their boundaries) may be found in polynomial 
time. Within the interior of each piece, G is either always 
true or always false. Then the disjunctive normal form is 
simply \J p{x £ int(P)}, where the disjunction is taken over 
all polyhedral pieces P (including lower-dimensional ones) 
such that G is true. □ 



The counting and generating function forms of this problem 
(either counting the number of solutions to or finding the 
generating function of a quantifier-free Presburger formula) 
were proven to be polynomial time by Barvinok [T] . 

Adding one quantifier alternation makes this problem hard, 
unfortunately: even with two variables, deciding 3aV& G{a, b) 
is NP-hard, as Schoning showed [55] • A possible modifica- 
tion is to also fix the number of linear inequalities allowed. 
In this case, Kannan proved [2Z! that the decision problem 
3aVb G(a, b) can be solved in polynomial time, and Barvi- 
nok and Woods proved ^ that the counting and generating 
function problems could be solved in polynomial time. 

With more quantifier alternation, the computational com- 
plexity of the problem is unknown. 



Conjecture 9. For a fixed number of variables and lin- 
ear inequalities, there exists a polynomial time algorithm to 
compute the generating function for a Presburger set. Hence 
the counting problem and the decision problem may be an- 
swered in polynomial time. 

5.2 B^c 

The first question here is how to interpret the problem: how 
should a quasi-polynomial be represented? The obvious way, 
listing every constituent polynomial, is inefficient. For ex- 
ample, the generating function 1/(1 — x") corresponds to 



the quasi-polynomial g{x) — 1 if njs and g{x) = oth- 
erwise; this requires n polynomials to represent. The so- 
lution of Verdoolaege and Woods [H] is to use piecewise 
step-polynomials: 



Definition 9. A step-polynomial g : N" — Q is a function 
written in the form 

m 

(p) = ^ Cti Y\ lO-ijlPl H ^ aijnPn + bij \ , 

,: = 1 j = l 

where Oi, aijk,bij G Q and [-J is the greatest integer func- 
tion. 



For example, the counting function corresponding to 1/(1 
x") can be written at the step-polynomial I — I — I ^ I • 



Definition 10. A piecewise step-polynomial is a function 
g : N" — > Q such that there exists a finite partition [J^iPi n 
N") of N" with Pi polyhedra (which may not all be fuU 
dimensional) and there exists step-polynomials gi such that 



5(P) = <7«(P) forpeP, nN". 



Verdoolaege and Woods prove (44) that one may convert 
between rational generating functions and piecewise step- 
polynomials in polynomial time (for fixed degree of the step- 
polynomial). Therefore, the conceptual equivalence, B ^ 
C, is also an algorithmic equivalence, provided piecewise 
quasi-polynomials are suitable represented. 

5.3 Other Complexity Relationships 

One reason generating functions are so valuable is that al- 
gebraic manipulations can be brought into play. We've al- 
ready seen that substituting in 1 for variables can be done in 
polynomial time, enabling us to compute cardinality. Sev- 
eral other interesting operations have polynomial time algo- 
rithms. 

For example, given the generating functions for two sets S 
and T, the generating functions for S CiT, S U T, and 5* \ 
T can all be computed in polynomial time 5 (assuming 
fixed number of variables and fixed bound on the number of 
binomials in the denominators of the generating functions). 
This means that the boolean operations and/or/not can be 
mirrored by generating functions. The key ingredient is the 
Hadamard product, which can be computed in polynomial 
time [5]. 



Definition 11. The Hadamard product of 

/(x) = c(p)xP and <?(x) = ^ d(p)xP 



pel* 



p6N" 



(/*g)(x)= Y. ^(pKpK 

pSN" 



For example, the fact that /(STlT; x) = f{S; x)*/(T; x) im- 
mediately shows us how to compute the generating function 
for SnT, given /(S';x) and /(T;x). 

Unfortunately, the operation of quantification cannot be sim- 
ilarly mirrored efficiently with generating functions: if we are 
given /(5';x, J/) for some S C N''^^, and define T = {a e 
N'' : V6 G N, (a, b) £ S}, then obtaining /(T; x) is NP-hard. 
To see this, take S to be the Presburger set for a quantifier 
free formula, G{a,b). We can compute f{S;x,y) in polyno- 
mial time, as discussed in Section [5.11 If we could compute 
f{T; x) in polynomial time, then we would be able to decide 
3aV6 G{a, b) in polynomial time, by checking if f{T; 1) 7^ 0. 
But deciding this is NP-hard l36l, as mentioned in Section 

Even simpler sounding tasks can be surprisingly difficult. 
There is no presently-known "nice" way to tell whether a ra- 
tional generating function, (/(x) (input as a sum of the basic 
rational functions from Definition [S]), is identically zero. A 
heavy-handed - but polynomial time - algorithm is, after 
checking that 1 is not a pole of g, to substitute 1 into g*g: 
the coefficients oi g * g are nonnegative, so (g* (;)(!) = if 
and only if g — 0. 
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